Submitted to: ACTA PHYSICA POLONICA A 



Magnetic orderings and phase separations in the zero-bandwidth limit 
of the extended Hubbard model with intersite magnetic interactions 
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A simple effective model for a description of magnetically ordered insulators is analysed. The tight 
binding Hamiltonian consists of the effective on-site interaction ([/) and intersite magnetic exchange 
interactions (J^, J^^) between nearest-neighbours. The phase diagrams of this model have been 
determined within the variational approach, which treats the on-site interaction term exactly and 
the intersite interactions within the mean-field approximation. We show that, depending on the 
values of interaction parameters and the electron concentration, the system can exhibit not only 
homogeneous phases: (anti-)ferromagnetic (Fa) and nonordered (NO), but also phase separated 
states (PSa: Fq-NO). 

PACS numbers: 71.10.Fd, 75.10.-b, 75.30.Gw, 64.75. Gh, 71.10.Hf 



I. INTRODUCTION 

The extended Hubbard model with anisotropic spin 
exchange interactions fl'-'S: is a conceptually simple phe- 
nomenological model for studying correlations and for 
a description of magnetism and other types of electron or- 
derings in narrow band systems with easy-plane or easy- 
axis magnetic anisotropy. 

In this report we will focus on the zero-bandwidth limit 
of the extended Hubbard model with magnetic interac- 
tions for the case of arbitrary electron density < n < 2. 
We consider the U-J^ Hamiltonian of the following form: 
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where U is the on-site density interaction, J^ is z- 
component of the intersite magnetic exchange interac- 
tion, ^/j ^-.j restricts the summation to nearest neigh- 



bours. 






denotes the creation operator of an electron 



with spin a at the site i, fii — '^^ fiia-, fii^ — cl^Cic and 
sf — \{hi'^ — fiii). The chemical potential n depending 
on the concentration of electrons is calculated from 
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with < 71 < 2 and N is the total number of lattice sites. 
The model (fTl) can be treated as an effective model 
of magnetically ordered insulators. The interactions U 
and J^ will be assumed to include all the possible contri- 
butions and renormalizations like those coming from the 
strong electron-phonon coupling or from the coupling be- 
tween electrons and other electronic subsystems in solid 
or chemical complexes. In such a general case arbitrary 
values and signs of U are important to consider. We re- 
strict ourselves to the case of positive J^ > 0, because of 



* |kakonrad@amu.edu.pl| 



the symmetry between ferromagnetic ( J^ > 0) and anti- 
ferromagnetic ( J^ < 0) case for lattice consisting of two 
interpenetrating sublattices such as for example sc or bcc 
lattices. 

We have performed extensive study of the phase dia- 
gram of the model ([T| for arbitrary n and jj, \W,T' . In the 
analysis we have adopted a variational approach (VA) 
which treats the on-site interaction U exactly and the 
intersite interaction J^ within the mean-field approxi- 
mation (MFA). We restrict ourselves to the case of the 
positive J^, as it was mentioned above. 

Let us point out that in the MFA, which does 
not take into account collective excitations, one ob- 
tains the same results for the [/-J^ model and the 
[/- J^^ model, where the term 2 J^ ^ sfs^ is replaced 

with J'^^J2i^tK^^tK)' describing interactions be- 
tween xy-components of spins at neighbouring sites, 
sf = ct-Cii = {s~)^. In both cases the self-consistent 
equations have the same form, only the replacement 
jz _j. jxy jg jjgg(jg(j and a magnetization along the z-axis 
becomes a magnetization in the xy-plane [6_ . 

For the model (111 only the ground state phase dia- 
gram as a function of yit [S] and special cases of half-filling 
(n — 1) [9J and f/ — >■ oo [TU] have been investigated till 
now. 

Within the VA the intersite interactions are decou- 
pled within the MFA, what let us find a free energy 
per site f{n). The condition (pi) for the electron con- 
centration and a minimization of f{n) with respect to 
the magnetic-order parameter lead to a set of two self- 
consistent equations (for homogeneous phases), which are 
solved numerically. The order parameter is defined as 
to" = (1/2)(to^ -I- m%), where m^ = ^ '£^e^ (s?) is the 
average magnetization in a sublattice j = A, B in the 
a = z,xy direction (s^^ corresponds sf here). If to," is 
non-zero the ferromagnetic phase (Fq,) is a solution, oth- 
erwise the non-ordered phase (NO) occurs. 

Phase separation (PS) is a state in which two domains 
with different electron concentration exist in the system 
(coexistence of two homogeneous phases) . The free ener- 



gies of the PS states are calculated from the expression: 

fps{n+,n_) =mf+{n+) + (1 -m)f_{n_), (3) 

where f± {n± ) are values of a free energy at n± corre- 
sponding to the lowest energy homogeneous solutions and 

_ "-"- jg ,^ fraction of the system with a charge den- 
sity «_)_. We find numerically the minimum of fps with 
respect to n+ and n_. 

In the model considered only PSq, state (i. e. a coexis- 
tence of Fq and NO phases) can occur. 

In the paper we have used the following conven- 
tion. A second (first) order transition is a transition 
between homogeneous phases with a (dis-) continuous 
change of the order parameter at the transition temper- 
ature. A transition between homogeneous phase and PS 
state is symbolically named as a "third order" transition. 
During this transition a size of one domain in the PS state 
decreases continuously to zero at the transition temper- 
ature. 

Second order transitions are denoted by solid lines on 
phase diagrams, dotted curves denote first order transi- 
tions and dashed lines correspond to the "third order" 
transitions. We also introduce the following denotation: 
Jq = ziJ" for a = z, xy, where zi is the number of near- 
est neighbours. 

Obtained phase diagrams are symmetric with respect 
to half-filling because of the particle-hole symmetry of 
the Hamiltonian (fl]), so the diagrams will be presented 
only in the range < n < 1. 



II. RESULTS AND DISCUSSION 
A. The ground state 

In the ground state the energies of homogeneous 
phases have the form: for NO: Enq — {l/2)Un 
and for F„: Ep = -{l/2)J^n^ if n < 1 and 
Ep = U{n - 1) - (1/2) Jo"(2 - n)2 if n > 1. Com- 

paring the energies we obtain diagram shown in Fig. IT] 
At U = ~Jq{1 — \n — 1|) the first order transition 
Fq-NO takes place in the system. This transition is 
associated with a discontinuous disappearance of the 
magnetization. 

The first derivative of the chemical potential 
d^/dn = d^E/dn? for U/Jq > — 1 in the lowest energy 
phases is negative what implies that homogeneous phases 
are not stable (except n = 1). 



B. Finite temperatures 

Finite temperature phase diagrams taking into account 
only homogeneous phases and plotted as a function of 
U/Jq for chosen n are shown in Fig. 12k. The tricriti- 
cal point Ti , which is connected with a change of tran- 
sition order, for n = 1 is located at ksT/J^ — 1/3 and 
U/Jg' = -2/3 In 2 m. 
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FIG. 1. Ground state phase diagrams as a function of n with- 
out consideration of PS states. The dotted line denotes dis- 
continuous transition. 



The range of the occurrence of Fq, phase is reduced 
with decreasing n. For n > 0.67 and any U/Jq > —1 
we observe only one transition F^-NO with increas- 
ing temperature. In the range 0.67 < n < 1 the U/Jq 
coordinate of the Ti-point remains constant, so for 
U/Jq < —2/3 In 2 the F^-NO transition is discontinuous. 
However, for n < 0.67 in some range of U/Jq there can 
appear a sequence of two transitions: NO-F^^-NO. 

In Fig. (2)3 there are shown dependencies of the tran- 
sition temperature F^j-NO as a function of n for chosen 
values of U/Jq. The range of F^ stability is reduced with 
decreasing of U/Jq. For U/Jq>0 and any n we observe 
only one second order transition F^-NO with increas- 
ing temperature. There exist ranges of n and U/Jq < 0, 
where the sequence of transitions: NO-Fq,-NO is present. 

At sufficiently low temperatures homogeneous phases 
are not states with the lowest free energy and there PS 
state can occur. On the phase diagrams, where we con- 
sidered the possibility of appearance of the PS states, 
there is a second order line at high temperatures, sep- 
arating Fq, and NO phases. A "third order" transition 
takes place at lower temperatures, leading to a PS into 
Fq and NO phases. The critical point for the phase sep- 
aration (denoted as T2, a tricritical point) lies on the 
second order line Fq,-NO. Phase diagrams for U/Jq = 1 
and U/Jq = 10 are shown in Fig. l3] 

In the ranges of PS stability the homogeneous 
phases can be metastable (if dji/dn > 0) or unstable (if 
d^/dn < 0). We leave a deeper analyses of meta- and 
unstable states to future publications. 
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FIG. 2. Phase diagrams (a) kBT/J" vs. U/Jq for fixed n and (b) ksT/jQ vs. n for fixed U/Jg without the consideration of 
PS states. Dotted and sofid fines denote first and second order transitions, respectively. 
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FIG. 3. Phase diagrams ksT/Jo vs. n with the consideration of PS states for: U/Jq = 1 (a) and U/J" = 10 (b). Sofid and 
dashed fines indicate second order and "third order" boundaries, respectively. 



III. FINAL REMARKS 

We considered a simple model for magnetically ordered 
insulators. It was shown that at the sufficiently low tem- 
peratures homogeneous phases do not exist and the states 
with phase separation are states with the lowest free en- 
ergy. On phase diagrams we also observe the tricritical 
points, which are associated with a change of transition 
order (Ti-point, Fig.l2| or are located in the place where 



the second order line connects with "third order" lines 
(T2-point, Fig.[|). 

Let us stress that the knowledge of the zero-bandwidth 
limit can be used as starting point for a perturbation ex- 
pansion in powers of the hopping and as an important 
test for various approximate approaches (like dynami- 
cal MFA) analyzing the corresponding finite bandwidth 
models. 
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